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Abstract. Within the framework of either Euclidian (functional- integral) quantum 
gravity or canonical general relativity the signature of the manifold is a priori 
unconstrained. Furthermore, recent developments in the emergent spacctimc 
programme have led to a physically feasible implementation of (analogue) signature 
change events. This suggests that it is time to revisit the sometimes controversial topic 
of signature change in general relativity. Specifically, we shall focus on the behaviour 
of a quantum field defined on a manifold containing regions of different signature. We 
emphasise that, regardless of the underlying classical theory, there are severe problems 
associated with any quantum field theory residing on a signature-changing background. 
(Such as the production of what is naively an infinite number of particles, with an 
infinite energy density.) We show how the problem of quantum fields exposed to 
finite regions of Euclidean-signature (Riemannian) geometry has similarities with the 
quantum barrier penetration problem. Finally we raise the question as to whether 
signature change transitions could be fully understood and dynamically generated 
within (modified) classical general relativity, or whether they require the knowledge of 
a theory of quantum gravity. 
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1. Introduction 

Signature change in general relativity is a subject that over the years has sometimes 
attracted considerable heated controversy. Kinematically, the key point is that the 
choice of signature, the pattern of the signs of the eigenvalues of the matrix of metric 
coefficients, is both a coordinate invariant, and is an independent choice that is made 
simply for observational physical reasons. The signature does not seem to be constrained 
at the level of dynamics [U Ej. Neither the Einstein-Hilbert action nor the Einstein 
equations are tied to a specific choice of signature. This feature of general relativity 
becomes most obvious in two places: in Euclidean quantum gravity (where the spacetime 
manifold is indeed formally continued to Euclidean signature) [3j EJ |5j [6], and in 
the canonical Hamiltonian formulation of (ordinary Einstein-Hilbert) classical general 
relativity (where the overall signature of the spacetime is determined by the behaviour 
of a non- dynamical part of the metric related to the "lapse function") [7]. Indeed, it is 
always possible to locally pick coordinates to foliate any arbitrary manifold in the form 

ds 2 = -B dt 2 + gij (dx i + JV* dt) (dx j + N j dt). (1) 

Here A" is the shift vector, and B is simply some real function of position. If B 
is positive, then B = N 2 , where N is the lapse function, and the manifold has 
Lorentzian signature. In contrast, if B is negative, then B = —N 2 , where N is now 
the "Euclidean lapse" function, and the manifold has Riemannian signature. (The 
coordinates themselves, t and x l , always remain real.) 

In analogue spacetimes (emergent spacetimes) the quantity B is proportional to 
the (speed) 2 of whatever signal one is interested in studying 013]. Specifically, in BECs 
this is the (speed) 2 of the phonons, which is proportional to the so-called "scattering 
length", which in turn can be externally controlled and driven negative by using an 
external magnetic field to tune the system through a Feshbach resonance [9j [TOj, [TT] . 
Thus analogue spacetimes can provide explicit physical models for implementing (some 
of the key features of) signature change in laboratory systems — a circumstance that 
strongly indicates that it may be worthwhile to re-assess the whole notion of signature 
change in general relativity itself. 

Specifically, how does a quantum field theory (QFT) residing on a spacetime 
undergoing signature change react to the imposition of a signature change event? A 
specific and relatively simple example of this phenomenon has been carefully studied 
by Dray et al. [2], who considered massless conformally coupled scalar fields in a 
two-dimensional spacetime undergoing a simultaneous universe-wide signature change 
process of the form 

(1 + 1) (0 + 2) (1 + 1). (2) 

Here the (a + b) notation denotes a time dimensions and b spatial dimensions in a region 
of constant signature. Their key result can be slightly rephrased as follows: In this 
particular process, the Bogoliubov coefficients governing particle production are given 
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(up to unimportant phases) by the exact expressions: 

a = cosh (c k At) ; (3 = i sinh (c k At) ; (3) 

where At is the total (Euclidean) proper time (At = j N dt) for which the universe 
exhibits Euclidean signature (B < 0), k is the wavenumber of the mode under 
consideration, cq is the physical speed of light, and the conformal flatness of arbitrary 
(1 + 1) spacetimes was an essential ingredient in their original calculation. Additional 
background material can be found in references [121 EES EH] • An immediate result of 
this calculation is that both the total number and total energy of the particles generated 
in such (1 + 1) — > (0 + 2) — > (1 + 1) signature change events is infinite 0: 

N # = j\f3\ 2 dk = j |sinh(c A;Ar)| 2 dA; = oo; (4) 

E = (j)c k\(3\ 2 dk = jc^k |sinh(c fcAr)| 2 dk = oo. (5) 

This feature of the calculation is disturbing, and one possible interpretation of these 
results is that placing a QFT on a manifold actually prevents signature change - 
this is similar to the situation regarding topology change where for some time it was 
similarly felt that infinite particle production by topology change events might suppress 
them completely [IB] . (Current understanding is more subtle, and progress regarding 
topology change has been made by enlarging the "universe of discourse" to contain more 
general classes of spacetime [IT].) Below we shall generalise and modify the Dray et 
al. result in several ways: 

• First, it is useful to verify that the basic result is not significantly modified by 
going to (d + 1) dimensions, and in particular the physically interesting (3 + 1) 
dimensions. Though generally one no longer has conformal flatness, this condition 
can for calculational purposes be replaced by a condition that the universe is not 
expanding during the signature change event. (That is, the signature change event 
takes place on a hypersurface of zero extrinsic curvature.) 

• Second, we shall soon see that the inclusion of particle masses, which is important 
for any realistic attempt at model building based on this cosmological particle 
production effect, also leads to a closed- form exact result 



a 



cosh < \/m 2 CQ + Cq k 2 At > ; (6) 



(3 = i sinh ^m 2 4 + c 2 k 2 Ar j ; (7) 

where again the universe is taken to be non-expanding during the signature change 
event. 

| We shall use N# whenever we are talking about total particle number to avoid any possibility of 
confusion with the lapse function N. 
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• Third, in situations where the universe is undergoing expansion (or contraction) 
during the signature change event, while one cannot generically write down exact 
expressions for the Bogoliubov coefficients, there are nevertheless reliable WKB 
estimates for these coefficients. Namely, up to irrelevant phases: 



where a(r) is the scale factor of the universe as a function of Euclidean proper time, 
and the integral runs only over the Euclidean region E. 

• Fourth, we shall also consider the (nowadays rather popular) possibility of an ultra- 
high-energy breakdown of Lorentz invariance [TSJ [TjJJ |2U] , to see what happens to 
the particle production spectrum in this situation. If the dispersion relation in the 
absence of cosmological expansion is some arbitrary function u = Q(k), then in the 
presence of expansion the dispersion relation becomes u = Q(k/a(r)), and we shall 
(up to irrelevant phases) derive the estimates: 



In all of these situations we still see infinite particle production, N# = oo, and infinite 
energy being dumped into the QFT modes, E = oo. Thus the infinities encountered in 
the Dray et al. calculation are not specific to their particular toy model, but are instead 
generic features of Lorentzian— >■ Euclidean— s-Lorentzian signature change events, with a 
behaviour quite different from ordinary cosmological particle production [21] . This is a 
circumstance that certainly calls out for some careful and delicate analysis. 

For instance, a more drastic choice based on an analogy with BEC physics that 
additionally breaks time reparameterisation invariance, can be used to successfully 
regularise total particle production, though total energy release is still formally 
infinite [H]. We emphasise that in the current article we are working strictly within 
the confines of standard general relativity, and are not adopting any "analogue model" 
point of view such as considered in our earlier article [11]. (For general background on 
these analogue model approaches see [HI El [22] . ) This approach changes the entire focus 
and intent of the investigation, potentially making it more interesting to the bulk of the 
relativity and cosmology communities. Staying within the confines of general relativity 
means we are not working with a bi-metric theory and technically simplifies calculations 
as we have the full coordinate invariance of standard general relativity available to us. 
Indeed, strategic use of the full coordinate invariance of standard general relativity (and 
in particular the use of full time reparameterization invariance) is essential to keeping 
the calculations tractable. 







(10) 
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In the discussion below we shall first look at the kinematics of signature changing 
spacetimes, carefully specifying the precise nature of the problem we wish to solve (in 
particular the relevant PDE and junction conditions). We then consider certain exact 
results we can derive for non-expanding but signature changing spacetimes, follow this 
up with a discussion of signature change in an expanding universe, and then turn to the 
question of a possible ultra-high-energy breaking of Lorentz invariance. Finally in the 
conclusions we shall turn to the question of what this all tells us about the possibility 
of signature change in the observable universe - - is the apparent pathology of this 
curved-space QFT calculation telling us that we need to step outside the framework of 
semi-classical gravity? Is this calculation telling us something about quantum gravity 
itself? 

2. Kinematics of signature change in (d + 1) dimensions 

2.1. Coordinate independent properties of signature change events 

The matrix of metric coefficients g a b is a real symmetric (d + 1) x (d + 1) matrix which 
has (d + 1) eigenvalues. These eigenvalues are of course coordinate dependent, but the 
pattern of signs of the eigenvalues is independent of one's choice of coordinates — this 
is a classic result most typically referred to as Sylvester's theorem (or Sylvester's law of 
inertia) [23] . The pattern of signs of the eigenvalues of the matrix of metric coefficients 
is referred to as the signature of the manifold. If one sign is different from all the 
others, this is taken to qualitatively indicate the "time" direction and the manifold is 
said to have Lorentzian signature. If all signs are the same the manifold is said to have 
Euclidean signature. (Lorentzian and Euclidean signature are the two most important 
cases.) If some eigenvalues are zero the manifold is said to have degenerate signature. 
(We shall take our manifolds to be almost everywhere non-degenerate. Degeneracies, 
if present, will be confined to isolated hyper-surfaces of signature change. Note that 
due to the coordinate invariance of general relativity, a coordinate reparameterization 
can be made to obtain a continuous signature change event.) The remaining cases are 
said to have Kleinian signature [21], corresponding for instance to more than one time 
and more than one space directions. It is important to note that signature can (in 
the first instance) be defined in a purely kinematical way before even introducing any 
notion of dynamics — of course if we then want to generalise the discussion to asking 
questions about a possible dynamical origin or explanation for signature change events, 
then consideration of the dynamics is unavoidable. 

2.2. Signature change events in FLRW spacetimes 

We are particularly interested in spatially flat (k = 0) FLRW-like spacetimes, where 
space and time are already decomposed in an appropriate manner, and thus the shift 
vectors are not necessary, A" = 0. Standard FLRW geometries exhibit one degree of 
freedom, the scale factor; we shall write the scale-factor-squared as A(t), and in principle 
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allow A(t) to change sign. In addition, we allow the lapse-squared to be an arbitrary 
function B(t), also capable of flipping sign 

ds 2 = -B(t)dt 2 + A(t)dx 2 . (12) 

For B(t) = 1 > and A(t) = a(t) 2 > we recover the standard representation for a 
spatially flat (k = 0) FLRW universe. More generally, B(t) and A(t) are in principle 
allowed to go through zero and to change sign. The relevant metric, inverse metric, and 
metric determinant are 

g ab = diag(-5, A,A,.. .); g ab = diag ( (13) 



g = -BA d - = VBA*. (14) 

The coordinate speed of light, which we shall denote by c, (as opposed to the physical 
speed of light, which we shall denote by Co), is defined by 

c 2 = B/A. (15) 

At a Lorentzian— )■ Euclidean or Euclidean— s-Lorentzian signature change event: 

• A(t) can in principle change sign (as a function of t), but is fixed (in the sense that 
it is independent of how one coordinatises t). 

• B(t) can change sign (as a function of t), and is variable under reparameterisation 
of t = fit): 

B{t) dt 2 = B{i) di 2 B = B (^) . (16) 

• Note however that at any particular point in spacetime on the manifold, the sign 
of B is fixed (independent of how one coordinatises t). That is, for any arbitrary 
foliation, the sign of B is fixed for all points in the same hypersurface. 

There are in principle (at least) 3 distinct types of signature change: 

- + + ...-> + + + ... (17) 

- + + ...-)• ... (18) 

- + + ...-> + --... (19) 



depending on whether B, A, or both change sign l§|. Only the first two cases are truly 
instances of Lorentzian— j-Euclidean signature change, the third case corresponds to a 
sign flip in an overall conformal factor but does not change that nature of the light cones 
(there is still one "special" direction which we can call time, and d "other" directions 
we can call space). The second case involves a Big Bang singularity as A(t) — > 0, (and 
so a(t) — > 0), which while interesting in its own right is not germane to the discussion 

§ We shall not, in the current article, consider more exotic forms of signature change, such as Kleinien 
signature [24] . where one deals with two or generally more timclikc directions (and the relevant physical 
questions become ones of barrier penetration and reflection). 
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at hand. So we shall only be interested in the case where B(t) alone changes sign, that 
is when 

- + + ...->• + + + ...->- + + ... (20) 

This will be our fundamental model of a Lorentzian— >Euclidean— s-Lorentzian signature- 
change event. Because of the coordinate invariance of general relativity one can 
parametrise the t coordinate t — > i so that quite generically and without any loss of 
generality: 

• We can choose the coordinate i so that B satisfies: 

A d /B = l/e = ±l. (21) 

That is, 

e = sign[£(i)] = sign[£(t)]. (22) 

Equivalent ly, 

* = = *(*)'. (23) 

• This particular choice is made to make the d'Alembertian operator look simple in 
these coordinates. Specifically 

A d+1 = -^=d a (V=g g ab d b ___) . (24) 



This becomes 



Aj +' = ^ (!-] +t- (25) 



We also have 



(§) = V^jB = yC¥J\B\- 



(26) 



Equivalent ly 

'di 



\B\A~ d/2 . (27) 



• For the coordinate speed of light we then have 

c 2 = B/A c 2 = B/A = eA d -\ (28) 

This completes our specification of the kinematics of the spacetime geometry we are 
interested in. 



2. 3. Differential equation and junction conditions 

Note that we can always use the coordinate freedom of general relativity to set B = ±A d 
and thus apparently make the signature change "discontinuous". This is however 
"merely" a choice of coordinates. A careful analysis (see reference p~|, and related 
discussions in Dray et al. [2]) agree on the following procedure: To see what the junction 
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conditions are it is sufficient to work at the classical level and consider for instance a 
massive minimally coupled classical scalar field governed by the Klein-Gordon equation 
1 



'-9 



d a (V=g 9 ab d b 



m 2 (f). 



(29) 



Now rewrite the Klein-Gordon equation in the explicit form 

n ab 



d a diag -y / A d /B, VBA d ~ 2 , VBA d ~ 



d b (f) = yfBM m 2 0, (30) 



(31) 



that is 

- d t i^A d /B d t (^j + V • (VBA d - 2 V0^ = ^BA d m 2 <j>. 

Fourier transforming in space (but not time) this becomes 

d t (^A d /B dtA = - ^BA d m 2 + V BA d ~ 2 k 2 0. (32) 

Now use the coordinate invariance to transform t — > i and substitute back into the 
ODE, then 



d. 



VTA™ m 2 + VeA 2d ~ 2 k 2 



(33) 



that is, (certainly everywhere apart from those "apparently discontinuous" points where 
e changes sign), 



Ve 2 A 2d m 2 + Ve 2 A 2d - 2 k 2 



(34) 



dt{d- t (t>) = 

We can rewrite this as 

d\<$> = -e A d [m 2 + k 2 /A] 0, (35) 

where at worst we will have to apply some junction conditions at the points where 
e changes sign. This should be compared with the Dray et al. calculation [2], which 
corresponds to the special case m — 0, d — 1. Physically we can reinterpret this (still 
classical) calculation as a specific parametric oscillator 

= -u 2 0, (36) 

with 

Q 2 = eA d [m 2 + k 2 /A] . (37) 
There is a discontinuous jump in the "potential" at the point where e changes sign: 

• The Lorentzian regions are "classically allowed" . 

• The Euclidean region is "classically forbidden". 

• If you think of u 2 as a "potential", then this potential has infinitely sharp walls 
where e changes sign, and u 2 — > ±|o)| 2 . 

• If one reinterprets this as a quantum mechanical barrier penetration problem, then 
in the Euclidean region one is attempting to tunnel through a barrier with exactly 
one half of the energy one would classically need to just skim over the top of the 
"jump" . 
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At the barrier walls you will have to apply the Dray et al. boundary conditions [2J: 

[0] = 0; [n a d a <P} = 0. (38) 

This second condition can be written in several equivalent forms. For instance, 



dr 



0. 



dr = |c| dt = y/\B\/A dt 



or, using 
we also have 

But since we are assuming A is continuous this is equivalent to asserting 
f~A* deb 

or finally 

r/ui 

0. 



dt 



(39) 



(40) 



(41) 



(42) 



(43) 



This last form of the junction condition is the one that is actually easiest to use. (And 
this is the form of the junction condition that one might have most easily derived 
by simply integrating the differential equation across the signature change event.) 
Summarizing, the differential equation we need to solve is 

<9f0= -e \u\ 2 0, (44) 

subject to the boundary conditions 

[0] = 0; [d- t <t>] = 0. (45) 

Though for ease of presentation we have so far phrased things directly in terms of the 
PDE which governs the field oscillations, and so have implicitly phrased the analysis in a 
first-quantized wavefunction formalism, the same results can be obtained in a number of 
different ways {e.g. we could rephrase everything in terms of continuity of the classical 
field and classical momentum density, or we could perform a second-quantized mode 
analysis). These junction conditions also carry through for curved-spacetime QFTs. 

3. Non-expanding universe 

Ultimately the differential equation we need to solve is 







-e \u\ 2 



(46) 
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where in the current situation u is piecewise constant. Let e change sign at i_ and i + . 
Then in the three regions of interest we have 

t < t- : i/j = exp(-z|w|t); (47) 

fe (*_,?+): tjj = Aexp(+\u)\i) + Bexp(-\u)\i); (48) 

?>i_: V = aexp(+i|w|t) + £exp(-i|w|i). (49) 

Applying the junction conditions at the first signature changing event we have 

exp(— «|a)|t_) = + Aexp(+\u\i_) + £>exp(— |a)|t_); (50) 

exp(— z|a)|t_) = zAexp(+|a;|t_) — iB exp(— \uj\i_); (51) 

and from again applying the junction conditions at the second signature changing event 
we obtain 

-Aexp(+|a)|£+) + B exp(— |u>|£ + ) = a exp(+i\u\i + ) + (3 exp(—i\u\i + ); (52) 
ylexp(+|o}|t + ) — _Bexp(— \u)\i + ) = ia exp(+i\oj\i + ) — i/3 exp(— (53) 

These are easily seen to yield 
1 + i 

A = — - — exp( — |w|t_) exp(+i|w|t_); (54) 
1 — % 

B = — - — exp(+|a;|t_) exp(+z|d)|t_); (55) 

leading to 

a = cosh{ \cj\ (£+—£_)]} exp{ —i\Co\ (t + — £_) }; (56) 
= i sinh{ \u>\ (i + - £_) ] } exp{ (t+ + t_) }. (57) 

Up to irrelevant phases we can cast this as 

a ~ cosh{ \u\ At}; (58) 
f3 ~ zsinh{ |w| At}. (59) 

But from the definition of |a)|, and inserting suitable factors of Cq, 

\u\ = A d / 2 yjm?4 + c% k 2 /A. (60) 

Remember that in the current calculation A is a constant, A d / 2 At = At, and we 
furthermore can without loss of generality set A — > 1 by rescaling the spatial coordinates. 
Thus 



| af| = cosh <. yj m 2 CQ + Cq k 2 At > ; (61) 

\(3\ = sinh | yjm 2 cl + cl k 2 At | . (62) 
From a second-quantized viewpoint, this implies particle production given by 

N k = sinh 2 j Jm 2 4 + c 2 k 2 At \ . (63) 



These formulae generalize the results of Dray et al. [2] , in that they are now applicable 
in any number of spatial dimensions, and also allow for particle mass. 
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4. WKB estimate 

If A(t) is not a constant, corresponding to an expanding (or contracting) spatially flat 
FLRW universe that undergoes a signature change event, then for "slowly expanding" 
universes one can use WKB theory to give a reliable estimate for the classical / quantum 
field excitations. For a minimally coupled scalar in (d + 1) dimensions undergoing a 
L — >■ E — > L signature change transition we found, in equation (1351) that the Klein- 
Gordon equation reduces to: 

<j> = -eA d [m 2 + k 2 /A] 0. (64) 

Here a choice of e = +1 corresponds to Lorentzian regions of spacetime, and in the 
intermediate Euclidean region, e = — 1. To estimate the particle production using WKB 
theory, we seek an approximate solution to the Klein-Gordon equation in each region, 
and connect the solutions at each hypersurface of signature change by applying the 
Dray et al. junction conditions. That is, the scalar field (f) and its conjugate momentum 
dt4> are continuous at a hypersurface of signature change — see (l38l) and ( l43l) . 

4-1. Approximate solutions in the three regions 

We choose the initial Lorentzian region to extend from (— oo,a). The Euclidean region 
exists from the first hypersurface of signature change at i — a, extending until the second 
hypersurface of signature change i = b. The final Lorentzian region occurs from (b, +oo). 
Note this problem is mathematically identical to that of flux hitting a steep-walled 
potential well. In this case the regions of Lorentzian signature would correspond to 
the two classically allowed regions, (—00, a) and (b, +00) while the classically forbidden 
region corresponds to the Euclidean region. The signature changing hypersurfaces play 
the same role as (sharp wall) "turning points" if we are looking at a barrier penetration 
problem. 

• In the initial Lorentzian region, we start by defining ojl = vO 2 , and consider the 
approximate solution 

exp(— % f?uj L di) , 
4>{t < a) = , (65) 

with conjugate momentum 

4>(i < a) — +i^fujL exp ( —i f u L di] — <p. (66) 

Throughout the article, we use the notation 7 = d?y = dj/di. 

• In the final Lorentzian region, we again proceed by defining uj l — y/^P, and the 
scalar field takes the (approximate) form 

exp(+i [* ULdt) exp(— i f* uj L dt) , 
Mt >b)=a PV p_ 1 + p K 1 , (67) 
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with conjugate momentum 

4>(t > b) = ia y/ujL,exp ^+i J ^ u^dfj — i(3 y/uLexp ^—i J ujl^\ — -— 4> ■ (68) 

• In the Euclidean region, which extends from the first signature change 
hypersurface at t — a, up to the second signature change hypersurface at i — b, we 
start by defining co E = y/—L0 2 . Approximate solutions for the scalar field can be written 

as 

exp ( + f * u E di ) exp ( - f * uj E di ) 
(f> = A '- + B '- , (69) 



with the "conjugate momentum" in the Euclidean region being given by 

4> = A ^[uj e exp ^+ J u E d^j - B yfw E exp J u E d^j ~ ■ (70) 

4-2. Junction conditions at the two allow ed-to- forbidden transitions: 

We now apply the Dray et al. junction conditions at the two signature change events. 
Note that at the two junctions we have 

Wl(gT) = W£;(a + ); w l (gT) = UJ E (a + ); (71) 
u E (b~) = cu L (b + ); u E (b~) = co L (a + ). (72) 
Before the first signature changing event, the scalar field takes the form 

4,(i<a) = eM ~ if jy Lit) , (73) 



while just after the hypersurface of signature change, but still within the Euclidean 
region, 



exp ( + f * u E dt ) exp ( - f * u E dt 
(a<i<b) = A ^ L + B ^ '- . (74) 



We connect the scalar field across the signature change event by employing the junction 
conditions we expounded on previously, these are, that both the scalar field and 
conjugate momentum are continuous: [0] = and [0] = 0. Continuity of the scalar 
field gives us 

l = A + B, (75) 
while from the requirement that the conjugate momentum is continuous at a, we obtain 
i = A — B. (76) 

Therefore 

A= 1 -^; B= 1 —^. (77) 
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<f>(t>a) 



1 



exp ( + f*u E dt) i 



x 



I exp 
— x 



(76 



whence, near the second signature changing event 



<Ki<b) 



l+i 



exp ( + uj e dt) exp (— u E dt 



UJ, 



1 - i 
+ x 



exp (- j b a uj e dt) exp [J* Ue dt) 



UJ, 



(79) 



Define 6 = + J o^e dt, so that just before the second signature changing event we have 



l+i 



exp 



Si u Edt 



-e 



exp ( uj E dt 



50) 



We now need to match this to the scalar field in the final Lorentzian region, that is 

exp(+i Curdt) exp(— i f, Ur,dt) 
(j){t>b)=a V }L ' + (3 lK Jb ' 



Again, we will utilize the junction conditions to do so. From continuity of the scalar 
field, [(p] = 0, we obtain 



-e 



a + f3. 



2 2 

Continuity of the conjugate momentum, [0] = 0, gives us 



l + i 



-e 



i(a — (3). 



2 2 
4-3. Bogoliubov coefficients: 

Solving these two simultaneous linear equations for a and /3 we find 

a = cosh G; /3 — i sinh 0. 
Earlier we defined 6 = + J^u^dt, and so 



52) 



53) 



a = cosh < + / ue dt 



(3 = i sinh 



oje dt 



54) 



55) 



Recalling that u E = A{t) d (m 2 + k 2 / A(t)) , we find that the cosmological particle 
production for an expanding or contracting spatially flat FLRW universe with a region 
of Euclidean signature is governed by the Bogoliubov coefficient 



(3 = i sinh 
which we can also rewrite as 
(3 = i sinh 



A{t) d/2 ^m 2 + k 2 /A(t) di 



I ^m 2 + k 2 /A(t) y/\B(t)\ dt 

J E 



(86) 



(87) 
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In terms of our original t coordinate, we can write this as 
= i sinh I ^m 2 + k 2 /A(t) y/\B(t)\ dt 

J E 

Now introducing the physical scale factor via A(t) — > a(t) 2 , and proper (Euclidean) 
via dr = y/\B(t)\ dt = y/\B(i)\ di, we have 

j3 = i sinh / \J m 2 + k 2 /a(r) 2 dr 
Je 

Finally, inserting appropriate factors of Co, 
j3 — i sinh 



2 4 + k 2 c 2 /a(r) 2 dr 



which is our promised result. Similarly 



a = cosh 



J m 2 cl + k 2 cl/a(r) 2 dr 
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(88) 
time 

(89) 
(90) 
(91) 



5. Consistency check: E — |V for a rectangular barrier 

As a simple check on the calculation, it is a classic and quite standard result that for 
quantum tunnelling through a rectangular barrier one has (see for example [251 P- 79], 
or PSH2Z]) 

T = ^ - E \ (92) 

AE(V - E) + V 2 smh 2 y2m{y - E) L\ ' 

where V is the height of the barrier, L is its width, and E is the incident energy. If we 
now take the special case E = |V we have 

1 1 

(93) 



(94) 



l + sinh 2 [ v /2m(\/- J E) L] cosh 2 [^j2m(V - E) L] ' 
from which, using the standard equivalencies 

\a\ 2 ^l/T and \(3\ 2 ^ (1 - T)/T = R/T, 
we see 

H 2 cosh 2 [v/2m(V - E) L] = cosh 2 [K L] O cosh 2 { \u\ At} . (95) 

That is 

|a| <-> cosh { \uj\ At} , \(3\ o sinh { \u\ At} , 

completely in agreement with our exact and WKB calculations. 



(96) 



Signature change events: A challenge for quantum gravity? 



16 



6. Quantum gravity (phenomenology) and signature change events? 

There are several reasons for wishing to look at issues related to "quantum gravity 
phenomenology", reasons such as possible Lorentz symmetry breaking, the limits of 
validity of the semi-classical gravity approximation, and the question of back-reaction. 
Regarding Lorentz symmetry breaking — one particularly important reason for looking 
at this topic is that it is known, in BEC-based analogue spacetimes, that the same 
physics that is responsible for Lorentz symmetry breaking is also responsible for 
regularising particle production, and in fact rendering the number of particles produced 
finite. It would be interesting to see if Lorentz symmetry breaking regularises the 
particle production in this more abstract context. (In fact it does not, and in the 
present context we shall soon see that regular isat ion of the particle production requires 
more than just generic Lorentz symmetry breaking, it requires a particular form of 
Lorentz symmetry breaking.) Beyond this, the semi-classical gravity approximation, 
and the question of back-reaction, are central to the question of whether it is possible to 
define a dynamics for signature change events — at the present level of analysis we have 
treated signature change kinematically, not yet making any direct attempt at assigning a 
suitable dynamical theory to this process. While dynamical theories of signature change 
(or signature selection) have been considered in the past (e.g., see [281 [291 ESS EI]), the 
(infinite) particle production we have encountered in this article very strongly suggests 
that any realistic dynamical theory must be very strongly dependent on back-reaction, 
and that failure to include the possible back-reaction due to the specific processes 
discussed in this article will render any possible dynamical system moot. 

6.1. (d+1) massive minimally coupled scalar field with Lorentz breaking 

There are numerous reasons, typically based on quantum gravity phenomenology, for 
thinking that it might be interesting to consider the possibility of an ultra-high-energy 
breakdown of Lorentz invariance [TBI EE ED] • Indeed, let us adopt a Lorentz-breaking 
formalism similar to that of Jacobson [32]. (See also [33].) Consider a PDE of the form 




where A^ + i is the spacetime D'Alembertian while is a purely spatial D'Alembertian. 
We shall again adopt definitions (TT3l) (TT4|) . The relevant D'Alembertian operators are 
now specified by 






(97) 




(98) 




(99) 



Then the relevant PDE reduces to the explicit form 




(100) 



Vrm f(-a^v 2 



)0 



V BA d m 2 0, 
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which we rewrite as 

- d t (y^A d /B dtA + V (VBA d - 2 V<f>\ - V BA d F(-A~ 1 \/ 2 )(f) = \TBA d m 2 0, (101) 
implying 

dt (y/A d /B dt^j = - \ r BA d {m 2 + F(A- l k 2 )} + y/BA d ~ 2 k 2 <f). (102) 

Now use exactly the same transformation, reparameterising t — > i , (and so setting 
B = e A d ) to make the LHS simple. One can just transcribe the previous discussion 
with minor modifications: 

d 2 t< p = -e A d [m 2 + F(k 2 /A) + k 2 /A] 0. (103) 

That is, we now have 

io 2 = eA d [m 2 + F(k 2 /A) + k 2 /A] , (104) 

where even in the Lorentzian signature region the presence of the F(k 2 /A) term indicates 
the presence of an explicit breaking of Lorentz invariance. Despite the presence of this 
Lorentz breaking term, the formulae for the Bogoliubov coefficients are only slightly 
altered. Indeed, approximate WKB solutions for the mode functions in the Euclidean 
regime are now 

± J ^m 2 + k 2 /A + F{k 2 /A) A d/2 di 



exp 



/ m 2 + k 2 /A + F(k 2 /A) A d l 2 



(105) 



leading, in a manner completely analogous to the preceding calculation, to the WKB 
estimates (up to irrelevant phases) for the Bogoliubov coefficients: 



a 



P 



coshjy ^/m 2 + k 2 /A + F(k 2 /A) A d/2 dt J ; (106) 

zsinhjy Vm 2 + k 2 /A + F(k 2 /A) A d/2 dt j . (107) 

It is easy to check that this reproduces all the other cases we have considered in 
appropriate limits. Explicitly introducing the standard scale factor A(t) — > a(t) 2 , 
utilizing the dispersion relation in the form 

u = fi(Jfe) = ^/m 2 + k 2 + F{k 2 ), (108) 

and working in terms of (Euclidean) proper time r, we finally obtain the promised result: 

a w cosh { / n(k/a(r)) dr I ; (109) 



' E 



J Q(k/a(r)) dr } . 



Signature change events: A challenge for quantum gravity? 



18 



Note that the adoption of this particular type of Lorentz symmetry breaking does not 
ameliorate the production of infinite numbers of particles, and infinite energy in the 
QFT, that we have encountered in previous calculations: 



jsmti z { I n(k/a(r)) dr } d d k; (111) 





E « <p u(k/ aoo ) smh z { J n(k/a(r)) dr } d d k. (112) 

These infinities can be tracked back to the fact that in the current class of models we 
have retained full coordinate invariance for reparameterisations of the time variable, so 
that on transforming t — >■ t the value for B was specifically chosen to make it equal to 
A d , and so effectively eliminate it from the PDE for the mode functions. By invoking 
the time reparameterisation invariance present in general relativity, and then suitably 
choosing the time coordinate, we have greatly simplified the calculation (rendering it 
tractable) at the cost of limiting the class of Lorentz violations one can consider. 

This is emphatically not what happens in analogue spacetimes, where the presence 
of a background laboratory time explicitly breaks time reparameterisation invariance 
and greatly influences the physics. Consider, for example, the Bogoliubov coefficients 
as presented in [TT]. There the authors present a real-life model, one that is capable 
of undergoing signature changes in the geometry emergent from an ultra-cold gas of 
weakly interacting Bosons, and where the Bogoliubov coefficients are of the form 



a 







coshjy yJ-B m 2 - e 2 qp k A — B k 2 /A A d ' 2 dtj ; (113) 
isinhjy yJ-B m 2 - e 2 qp k* - B k 2 /A A d/2 dtj , (114) 



such that the extreme ultraviolet modes are unaffected by the sign change of B. (The 
constant e qp = h/(2M) in front of the k A term is in this situation determined by 
the physical mass M of the bosons that undergo condensation, and is emphatically 
not affected by signature change.) Sufficiently high energy modes will maintain their 
oscillatory behaviour, beyond the region of Riemannian signature. In some sense the 
introduction of Lorentz symmetry breaking in the current formalism has made the 
situation worse — to regularise the number and total energy of the particles produced 
in the signature change it seems that one must do something much more drastic than 
just break Lorentz symmetry: It seems that one needs to break time reparameterisation 
invariance, since only then can one hope to give an explicit and invariant meaning to the 
idea that the "amount" of signature change, quantified by B(t), might be "small". As 
long as one retains time reparameterisation invariance then the magnitude of B(t), 
though not its sign, is coordinate dependent; and asking that signature change be 
"small" is a meaningless request. 
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6.2. Dynamics of signature change events 

Regarding the possible dynamics of signature change events, it is important to realise 
that in standard general relativity the signature is not a local dynamical variable 

- at best it is a "non-local" and discrete dynamical variable associated with entire 
regions, rather than a dynamical field associated with specific points on the manifold. 
Specifically, in the standard Einstein-Hilbert implementation of general relativity the 
lapse and shift show up as Lagrange multipliers enforcing the super-Hamiltonian and 
super-momentum constraints, and so are explicitly non-dynamical. This observation 
may however be altered in certain modified theories of gravity, specifically in theories 
that break time reparameterisation invariance. Models such as BEC-based analogue 
spacetimes, because they are bi-metric, (they contain both a background spacetime 
metric and an "acoustic metric") explicitly break time reparameterisation invariance 

- the ratio between the speed of sound and the speed of light is a dimensionless 
observable that can be given an explicit coordinate independent meaning and this 
dimensionless ratio certainly has dynamics (dependent on the details of the BEC system 
under consideration). Effectively the lapse function of the acoustic metric has acquired 
a dynamics due to the presence of the background (laboratory) spacetime. 

If (and this is by no means certain) the true physics of "quantum gravity" is 
ultimately based on a bi-metric theory, then similarly the low-energy lapse function 
associated with "our" spacetime metric is likely to acquire dynamics, but only in those 
situations where the coupling of "our" spacetime metric to the "other" (background) 
spacetime metric becomes non-negligible. However — and this is important — given that 
our naive calculations (which ignored back-reaction) have resulted in infinite particle 
production and infinite energy release, it becomes clear that (barring miracles) any 
attempt at including back-reaction should lead to a strong coupling between the two 
metrics of any bi-metric theory — at which stage the signature becomes dynamical, one 
leaves the semi-classical regime, and the theory of "quantum gravity" must be invoked. 

It is in this sense that signature change is ultimately a challenge for the as yet 
unknown theory of quantum gravity itself - - attempting to analyse it within the 
framework of curved-spacetime QFT, or even semi-classical gravity, has led to an 
impasse. Like the questions related to topology change [IE] , and "chronology protection" 
(the avoidance of closed timelike curves) |34j . it seems that detailed analysis of this 
issue inexorably moves one outside the framework of semi-classical gravity it is 
ultimately "quantum gravity" (whatever that may be) that one will have to face in 
order to definitively answer the questions we raise in this article. 

7. Discussion and Conclusions 

The possibility of signature change in ordinary general relativity is a topic that in the 
past has lead to considerable interest jH EJ EJ EJ [121 CE] an d often controversy [H] C5] . In 
this article we have made a contribution to the theory of cosmological particle production 



Signature change events: A challenge for quantum gravity? 



20 



due to a Lorentzian— ^-Euclidean— >-Lorentzian signature change. 

• One of the key technical points of the article is the realization that cosmological 
particle production in signature change is closely related to quantum tunnelling 
"half-way up the barrier". 

• The WKB estimates we have derived can very easily be generalized to provide 
estimates for quantum tunneling through sharp walled barriers (with generic 
complicated peaks and valleys, as long as the allowed <H- forbidden transition occurs 
at the sharp vertical wall). 

• We have seen that generically signature change seems to drive the production of an 
infinite number of particles, with infinite total energy, and that these infinities are 
not ameliorated by dimension, rest mass, or even the most physically reasonable 
sub-class of Lorentz symmetry violations. 

• If one wishes to retain signature change as a physically viable process then it 
really looks like one should be searching for a suitable dynamics to both drive 
the signature change event and to regularise both the number and energy of the 
resulting particle emission. The only currently known mechanism for doing so 
is based on the presence of a preferred background time, of the type occurring 
in the so-called "analogue spacetimes". Of course the implied breakdown of 
time reparameterisation invariance will be deeply disturbing to classical general 
relativists. 

• One way of phrasing our results is this: Quantum field theory reacts violently to 
any attempt at imposing signature change on the spacetime manifold. In a manner 
reminiscent of the way that quantum field theory reacts violently to any attempt at 
imposing topology change [16], or chronology violations [M], imposing an externally 
driven signature change seems to catapult one into the realm of "quantum gravity" . 

Overall, we have seen that computational and conceptual advances in signature changing 
spacetimes can still be made, with what is from a technical perspective a remarkably 
straightforward and "fundamental" calculation. 
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